Excited-State Density Distributions of Neutron-Rich Unstable Nuclei 
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We calculate densities of excited states in the quasiparticle random-phase approximation (QRPA) 
with Skyrme interactions and volume pairing. We focus on low-energy peaks/bumps in the strength 
functions of a range of Ca, Ni, and Sn isotopes for = + , 1~, and 2 + . We define an "emitted- 
neutron number" , which we then use to distinguish localized states from scattering-like states. The 
degree of derealization either increases as the neutron-drip line is approached or stays high between 
the stability line and the drip line. In the 2 + channel of Sn, however, the low-lying states, not even 
counting surface vibrations, are still fairly well localized on average, even at the neutron drip line. 

PACS numbers: 21.10.Pc, 21.60.Jz 



The structure of excited states in exotic nuclei has been 
much studied recently, both in nuclei between Li and O 
PI 0j OJ 13 and in heavier Ca and Sn isotopes. Low- 
energy strength is often enhanced, and theorists have 
tried to understand the mechanisms responsible. Besides 
strength functions, they have examined transition densi- 
ties, which can be measured and provide information on 
excited-state structure. Our earlier paper contains an 
investigation of strength functions and transition densi- 
ties in a large range of medium-heavy spherical nuclei, 
from one drip line to the other. 

Transition densities tell us about where in the nucleus 
transitions occur; they reflect the spatial distribution of 
products of single-particle and single-hole wave functions 
and, as a consequence, are always localized, though they 
can be very extended near the neutron drip line. In nu- 
clei near stability large transition strength implies collec- 
tivity, which in turn tends to cause localization (as we 
shall see in giant resonances). Near the drip line, how- 
ever, the familiar connection between strength and col- 
lectivity is less systematic at low energies. Reference 
shows many examples of behavior like that characteriz- 
ing light halo nuclei and predicted for certain + states in 
heavy nuclei in Ref . : large strength coming from non- 
collective transitions to very spatially-extended single- 
particle states. Those excited states are usually unbound, 
and may not even be quasibound, that is they may be 
completely delocalized scattering states 0. 

In this letter we suggest a measure of localization and 
then investigate the extent to which the strong low-lying 
excited states are localized, independent of their collec- 
tivity. To uncover changes in structure near the drip 
line, we track the degree of localization as N increases 
in J T = + ,l~, and 2 + channels. Localization, which 
we analyze by examining diagonal density distributions 
of the excited states, has not been systematically studied 
before and offers a new window into excitations in exotic 
neutron-rich nuclei. 

We can adopt ideas from simple one-particle quantum 
mechanics to distinguish localized states from extended 



ones. If the energy of a single particle in, e.g., a square- 
well potential is negative (Chap. 3 of Ref. then the 
tail of its wave function decays exponentially, and if the 
energy is positive the tail oscillates. In certain small 
positive-energy windows, however, the amplitude of the 
oscillating tail is much smaller than that of the wave 
function inside the potential. Such states are sometimes 
called "quasibound", and the solutions that exhibit no 
enhancement inside are sometimes referred to as "scat- 
tering states." 

Of course we are dealing with many-body quantum 
mechanics and the asymptotic wave function depends on 
many coordinates. The tail of the one-body density, how- 
ever, should still give us a measure of localization. In this 
letter, we use the Skyrme-QRPA wave functions (with 
the parameter set SkM* and volume pairing) from Ref. 
0, in which we discussed transition densities, to obtain 
the diagonal density of excited-state k: 
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where q is proton or neutron, r{ is the position operator 
of particle i, ipK(f) is a single-particle wave function in 
the canonical basis, and uk and vk are the associated oc- 
cupation amplitudes. The QRPA amplitudes X^ K , and 
F| K , are assumed real, K denotes the state conjugate to 
K, and we have used the convention V^( r ) = ( r )- 
The summation is taken only for protons or neutrons, and 
in our numerical calculations we will use the Hartree-Fock 
(-Bogoliubov) approximation to the ground-state density 
Pg(r). To obtain Eq. JIJ we have used the excited-state 
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where a* K and ajf are the creation and annihilation op- 
erators for quasiparticles in the canonical state K. 

If, as in all the calculations presented here, the ground 
state is spherically symmetric, the number of particles be- 
tween two radii is just the integral over that well-defined 
region of p q k {r)r 2 . We define the "emitted-neutron num- 
ber" associated with the state k as 
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where r c denotes the radius at which the density pjj(r) 
begins to develop a scattering tail. If N e (k) is close to 
or larger than 1, then the excited state is scattering- like 
(not localized), and if N e (k) is appreciably smaller than 
1, then the excited state qualifies as quasibound (local- 
ized). Intermediate cases are also possible, of course. If 
the pairing gap is not zero p^(r) does not integrate to the 
correct particle number 1 ; we therefore multiply p^(r) by 
a constant (at most a few percent from unity) to normal- 
ize it correctly. 

Rather than focus on individual states, we want to 
examine peaks in the strength function, which can en- 
compass several discrete states. We therefore define a 
strength-weighted average of the emitted-neutron num- 
ber for the strongly excited states k within a given peak: 
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where Sj? and Sjy denote isoscalar and isovector tran- 
sition strengths to the state k. The number of terms in 
the sums is between 1 and 10. We average the isoscalar 
and isovector quantities because both channels appear in 
low-energy strength- function peaks, particularly near the 
drip line. The giant resonances do not have this prop- 
erty, however, and for them we will use Nl s or as 
measures. 

Before investigating the dependence of emitted- 
neutron number on N and Z , we need to see whether we 
can calculate it reliably. The tail of the nuclear density is 



1 The QRPA does not guarantee the conservation of a quantity 
that is of second order or higher in the X's and Y's. 
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FIG. 1: Strength distributions in the isoscalar + channel of 
50 Ca, with a box size 20 fm. The neutron chemical potential 
is -6.61 MeV. 



extraordinarily sensitive, and converges slowly both with 
the size of our spatial "box" and the number of canoni- 
cal states in our basis. Figure Q shows the isoscalar + 
strength distribution for 50 Ca in a few versions of our 
calculation. We have varied two cutoff parameters: an 
upper limit e cr it on canonical single-particle energies (for 
protons, which have a vanishing paring gap), and a lower 
limit w cr it on occupation amplitudes (for neutrons, which 
have a non- vanishing gap). We confirmed that when 
e cr it = 250 MeV and t> cr it = 10~ 14 , the strength distribu- 
tion is almost identical to that with e cr j t = 200 MeV and 
^crit = 10~ n , shown in the bottom panel. In Refs. 0,^3 
we folded the strength to account for the 20-fm box ra- 
dius, and the values £ crr t = 150 MeV and v a a = 10~ 6 
were sufficient for convergence of the strength function 
and transition densities to the level of accuracy we re- 
quired. Figure ^ seems to suggests, however, that larger 
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FIG. 2: Neutron-density distribution of the low-energy state 
with non- negligible strength in 50 Ca (see Fig. 0, calculated 
with e cr it = 150 MeV, ii C rit = 1CP 6 , and a box radius of 20 
fm. The excitation energy is at the top of the figure. The 
ground-state density distribution (gs) is also drawn. In this 
case, r c — 7 fm. 



spaces are required to reproduce the tiny part of the 
excited-state density that lies significantly outside the 
nucleus. 

To examine the issue more closely, we show in Fig. |2 
the neutron density distribution associated with a low- 
lying state in the top part of Fig. ^ that distribution 
was obtained with the cutoff parameters we used exten- 
sively in Ref. In Fig. [3] we show the neutron densities 
associated with the corresponding states in the bottom 
part of Fig. obtained in a larger space. The tails of 
the more accurate densities are noticeably different (and 
more realistic). Yet improving the description this way 
has a relatively small effect on the emitted-neutron num- 
ber. In the smaller space the N e of low-energy states 
is 0.84, with individual N e (k) ranging from 0.77 to 0.90, 
and in the larger space it is 0.89, with N e (k) ranging from 
0.85 to 0.92. These ranges, 10% or less of the average, are 
not only small, but also typical of cases with N e near 1. 
Changing the box radius to 25 fm fragments the strength 
by adding new states, so that the N e (k) cannot really be 
tracked as the box size changes. The average N e over the 
states in a certain energy range, however, turns out to be 
a fairly robust measure of localization; with a box size of 
25 fm, N e — 0.86 for 5 sample states. Perhaps the sta- 
bility of N e , an integrated quantity, is not so surprising 
given that the smaller space is large enough to describe 
the wave functions extremely well for r less than r c . In 
any event, we can conclude here — for the + bump be- 
low the giant resonance in 50 Ca — that the states are 
scattering-like. Nearly a full neutron is far outside the 
nucleus, no matter how we vary parameters. 

To set a benchmark of sorts, we have looked at densi- 
ties for giant resonances in the isoscalar and isovector + , 
l - , and 2 + channels of several Ca, Ni, and Sn isotopes 
near stability. N~l s and iVg V are always on the order of 
0.5 or less. The isoscalar giant dipole resonances have 
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FIG. 3: The same as Fig. H but for e crit = 200 MeV and 

Vcrit = 10~ n . 



the largest values. The smallest belongs to the isovec- 
tor giant quadrupole resonance of 58 Ni, with iV* = 0.01. 
The reasons for this very small value are that the main 
components of the neutron excitation involve deep hole 
states and particle states near the Fermi surface, and that 
protons contribute appreciably. 

Having established a good measure of localization and 
calibrated it for typical collective resonances, we are now 
in a position to examine the behavior of low-lying exci- 
tations in neutron-rich nuclei. In each channel we look 
at low-energy bumps in a range of Ca, Ni, and Sn iso- 
topes. The resulting N e , plotted versus N, appear in 
Fig. QJ For + states they are more or less constant, at 
least in Ca and Sn, and since the constant is greater than 
0.8, the wave functions in the bumps are not localized. 
The N e in the 1~ channel in Ca and Ni show a clear 
increase with N , from 0.1 to about 1; those states are 
localized near stability and become more scattering-like 
towards the drip line. The 1~ states in Sn are harder to 
interpret; the curve has a maximum of about 0.8 around 
N = 100, but shows no clear trend in one direction or the 
other. The 2 + states in Ca and Ni mirror the behavior of 
the 1~ states. In the 2 + channel of Sn, N e increases with 
N but the low-lying states never get less localized than 
typical giant resonances. By our criterion, those states 
are localized, on average, even at the drip line, though 
they fluctuate within a given peak in a way that can- 
not be seen in the average (some individual states have 
N e (k) > 0.8). We should note that the states used to de- 
termine N e do not include low-lying surface quadrupole 
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vibrations, which usually have N e (k) < 0.1. 



1.6 
1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 


1.6 
1.4 
1.2 

1 

0.8 
0.6 
0.4 
0.2 


1.6 
1.4 
1.2 

1 

0.8 
0.6 
0.4 
0.2 





I 1 I I 1 1 I 1 I 1 1 I 1 1 I 1 I I 1 I 1 1 I 

: o + 


i i i | 1 1 i i i i i i i | i i i i i i i i i | i i i i 

Ca ■ : 

Ni — ^— - - 
Sn • . 


■_ h <3 

1 1 


• • 

■ ■ 1 1 I 1 1 1 1 i 


: r 




- , p 

- o< , 

1 1 


- 
'•- 

•. 

~* 

■ ■ 1 1 I 1 1 1 1 i 


: 2 + 


■ 






/ ? 


• 


'- / -• a ' 

Ml i ■ ■ 1 ■ i ■ i^i ■ i ■ 1 


'■ 1 i 1 1 ■ 1 1 


40 60 


80 100 120 



In Ref. |?J we found that the strength to low-lying 
bumps grew wth N in essentially every channel and every 
isotopic chain. Though we sometimes see similar behav- 
ior in N e , it is not universal. The + states, for instance, 
are not localized even near stability. And the 2 + states 
in Sn, while they become more scattering-like with in- 
creasing N, do not always delocalize beyond the level of 
giant resonances. All this comes with a caveat, however: 
it is possible that many-particle many-hole states that 
admix with QRPA excitations can alter N e somewhat. 
The one-particle one-hole nature of our states means, for 
instance, that N e is never much greater than one. Cal- 
culating the structure of many-particle many-hole states 
is difficult, however, and it may be some time before a 
systematic study of the kind reported here can be made. 

To summarize quickly: giant resonances consistently 
have N e around 0.5. In low-energy peaks N e is more 
varied, and changes with neutron number in ways that 
depend both on the isotopic chain and on multipolarity. 
Stronlgy excited low-energy states near the drip line are 
sometimes scattering states, but not always. 
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FIG. 4: N e in low-lying strength-function peaks, versus N, 
for the + , 1~, and 2 + channels of a sequence of Ca, Ni, 
and Sn isotopes. In the 2 + channel, the lowest-energy states 
(surface vibrations), carrying the largest strengths, are not 
included. In Ni the 1~ curve has more points than those 
of the other J 7 ' channels because the low-energy bump first 
appears at smaller N. 
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